Universality of the negativity in the Lipkin-Meshkov-Glick model 
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The entanglement between noncomplementary blocks of a many-body system, where a part of the 
system forms an ignored environment, is a largely untouched problem without analytic results. We 
rectify this gap by studying the logarithmic negativity between two macroscopic sets of spins in an 
arbitrary tripartition of a collection of mutually interacting spins described by the Lipkin-Meshkov- 
Glick Hamiltonian. This entanglement measure is found to be finite and universal at the critical 
point for any tripartition whereas it diverges for a bipartition. In this limiting case, we show that 
it behaves as the entanglement entropy, suggesting a deep relation between the scaling exponents 
of these two independently defined quantities which may be valid for other systems. 
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I. INTRODUCTION 

Considerable attention has recently been devoted to 
the study of genuinely "quantum" correlations or en- 
tanglement in the ground states of many-body systems 
among theorists [1] and experimentalists [2], an exciting 
field that profits from the interplay of quantum infor- 
mation and condensed matter. The measure of entangle- 
ment most extensively studied so far is the von Neumann 
entropy £ , which quantifies the entanglement between 
two complementary parts of a system, the common state 
of which is pure. This measure is known to display a 
universal divergence at criticality [1]. However, the sce- 
nario of "complementary parts" is rather restrictive by 
demanding that the two parts span the whole of a many- 
body system. Much more natural is the case where a 
portion of the system does not belong to either of the 
parts under consideration; that is, it forms an environ- 
ment. Moreover, such a general setting is of relevance 
for a situation where two parties have access to limited 
groups of constituents and attempt to exploit the entan- 
glement between these groups for quantum information 
tasks. 

In general, the state of two noncomplementary parts 
p is a statistical mixture and £ is no longer suitable to 
quantify their entanglement. To this end, one must in- 
voke the logarithmic negativity C [3], which is the only 
adequate measure for this task which is, at the same time, 
straightforwardly computable. In addition, it has oper- 
ational meaning (in terms of bounds) in teleportation 
and distillation. It is defined as C = lnTr((p Tl )' l '/5 Tl ) 1 / 2 , 
where T\ denotes the partial transposition that amounts 
to (m,n\'p Tl \k,l) = (k, n\p\m, I) with respect to a com- 
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plete set of basis states \k,l) = \k) 1 <E)\l) 2 on the bipartite 
system. Though, for the special case of pure states C and 
£ can both be computed from the so-called Schmidt coef- 
ficients [3], whether the former is universal and whether 
its behavior can be related to scaling exponents are open 
questions. Indeed C is defined purely from quantum in- 
formation considerations, and its manifestation of the el- 
egant scaling features from many-body physics will be a 
true surprise. 

Unfortunately, it is notoriously difficult to compute C 
even for the simplest of one-dimensional (ID) models for 
which some numerical results have been obtained [4, 5]. 
In this article, we present an analytic study of entangle- 
ment, as measured by C, between two macroscopic groups 
of spins, in a tripartite splitting of a many-body system 
displaying a quantum phase transition. We mainly focus 
on the Lipkin-Meshkov-Glick (LMG) model for which we 
find that C is always finite in stark contrast with other 
macroscopic correlation measures such as mutual infor- 
mation T which diverges. Remarkably, C does not de- 
pend on the anisotropy parameter at the transition point 
and may, in this sense, be considered as universal. How- 
ever, we also found that this is not a generic feature since 
the same analysis in the Dicke model [6] shows a different 
behavior. Most importantly, we show that in the limit- 
ing case of a bipartition, C diverges at the critical point 
as 1/6 In N ( where N is the system size ), exactly as the 
entanglement entropy. This property which is also found 
in other models leads us to conjecture that it should be 
valid for all systems. 



II. THE MODEL 

Let us consider a system of N spins 1/2 which are 
mutually coupled through an anisotropic XK-type inter- 
action and subjected to a magnetic field of strength h 
pointing in the z direction. The ground-state entangle- 
ment of this model introduced by Lipkin, Meshkov, and 
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Glick (LMG) in 1965 [7-9] to describe nuclei, has been 
the subject of many recent studies [10-18]. The LMG 
Hamiltonian is given by 



H=--(S 2 x + 7 S 2 y )-hS z , 



where S a = X^feLi <J 'a > /"^ ( a — x ,Ui z )i Co?' being the 
Pauli operators acting on the state space of the fcth spin. 
Here, we only consider the case of ferromagnetic inter- 
actions and, without loss of generality, we restrict, in a 
first step, the anisotropy parameter to ^ 7 < 1 and 
the field to h 0. This system undergoes a second-order 
quantum phase transition at h = 1, between a symmet- 
ric (h > 1) and a broken (h < 1) phase, which is well 
described by a mean-field approach. The corresponding 
classical ground state is fully polarized in the field direc- 
tion ((cr*) = 1) for h > 1, and twofold degenerate with 
(a l z ) = h for h< 1 (see Refs. [12, 19] for details). 
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III. £ IN A TRIPARTITION 

In order to compute entanglement of the ground state 
in a tripartite setting, we divide the TV spins into three 
groups 1, 2, and 3 with N±, N2, and N 3 spins, respec- 
tively, satisfying N1+N2+N3 = N. Accordingly, we par- 



tition the spin operators into S a = Sa^ + Sa' + Sa' ■ To 
diagonalize the Hamiltonian and hence obtain the ground 
state in the thermodynamical limit, it is convenient to ex- 

(k) 

press the S a 's in terms of bosonic operators using the 
so-called Holstcin-Primakoff representation [20] 



?(2) 



: (3) 



S<*> = N k /2 



a[a k , 



a\a k /N k 



,1/2 



(Sf )t, 



(2) 
(3) 



where = si^ ± iSy k \ Note that we focus here 

on maximum spin sectors to which the ground state of 
Hamiltonian (1) is known to belong. In addition, as dis- 
cussed in Ref. [12] for the single-mode case, substituting 
the expressions (2) and (3) in (1) requires a prior rotation 
of spin operators to bring the z axis along the classical 
magnetization of the ground state. Using this bosonic 
representation and expanding H at order (1/N k )°, one 
obtains 
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= 4 A k,l ai + l (4 B k ,i 4 + H.c.) + Ctc, (4) 

fc.Z=l 



where Cte denotes constant terms which will be irrele- 
vant in the following and where we introduced the 3x3 
real symmetric matrices 

Tl \fT\T-2 y/TlT 3 \ 

I Vn^2 t 2 jWfi , (5) 

y/TlT 3 y/T 2 T 3 T 3 I 



with Tfc = N k /N and I denotes the identity matrix. In 
the symmetric (h ^ 1) and broken (0 ^ h < 1) phases 
the prefactors read 



7-1 

-f — h 2 



0^h<l, 
€ h< 1 , 



(6) 



(7) 



and we further note that r > and r > s. Since the 
Hamiltonian (4) is quadratic, it is straightforwardly di- 
agonalizcd via a Bogoliubov transform. As discussed in 
[11] within a single-mode description, the gap for h ^ 1 
is given by A = \/(h— l)(h — 7), whereas it vanishes as 
cxp (— N) in the broken phase [19]. 

Without loss of generality, we compute the entan- 
glement of the ground state between group 1 and 
group 3 by computing the logarithmic negativity C be- 
tween the corresponding bosonic modes. To capture the 
entanglement properties of the arising mixed state the 
only available measure is C except for the limit of two 
spins for which concurrence can also be used [1]. In 
the present context a convenient definition for £ may 
be invoked. Indeed, the ground state is a Gaussian 
state for which an elegant framework is available that 
has allowed for a multitude of significant analytical re- 
sults in the past (see, e.g., [21]). Entanglement between 
any two of the groups may be inferred from the co- 
variance matrix T which collects the second moments 
T id = (il> \RiRj +RjRi\^o), (withi,j = 1,2,..., 6) of 
canonical coordinates x k = (4 + a k)/V^ and momenta 
Pk = i(4 — a k)/V%, which we group together in the 
vector R = (xx, x 2 , x 3 ,pi,p 2 ,p 3 ). In this representa- 
tion r adopts an explicit expression in terms of ma- 
trices V x = (A + B) and V p = (A — B) [22] which, 

by virtue of [14,^] = 0, reads T = r x ® r p , where 

_ 1 1 

T x = Tp 1 = V p 2 Vx . Similarly, the reduced density op- 
erator p = Tr 2 (|^o)(^o|) (Tr 2 denotes the partial trace 
over group 2) has a representation in terms of a Gaus- 
sian state with covariance matrix T obtained from Y upon 
canceling rows and columns that correspond to mode 2, 
namely, F = T x © T p , where 



T x =I+(a- 1 -1) 



n v /T i T 3 

^Jt\t 3 t 3 



r„=I+(a-l)( 



(8) 
(9) 



and a = y/(r + s)/(r — s) > 0. Then, upon partial 
transposition T\, the covariance matrix is subjected to 
partial time reversal p\ — > — <p\ [23] and is transformed 
into r Tl obtained from T by changing the sign in the 
off-diagonal terms in Eq. (9). The logarithmic negativity 
can then be computed in terms of the symplectic (degen- 
erate) eigenvalues of T Tl , 



Al,: 



1 + 9 ± y/g 2 +4riT 3 (a + a- 1 - 2), (10) 
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where g = [n + t 3 - (n - r 3 ) 2 ] (a + a -1 - 2) /2 > and 
where the + (— ) sign corresponds to Ai (A 2 ). Indeed, 
noting that Ai ^ 1 > A 2 > 0, the logarithmic negativity 
reads as [3] 



lnmin(Ai, 1) 



In A? 



(11) 



Equation (11) is the central result of this article and we 
now discuss it in detail. First, let us note that for an 
arbitrary tripartition of the system (0 < Tfe < 1, for 
k = 1,2,3) C is finite across the whole phase diagram 
(/i,7) including the transition point h = 1, where one 
has 



£(n,r 3 ,7, h = 1) 



1, n + r 3 - (n + r 3 ) 2 
- In — 
2 n 



73 - (ri - r 3 ) 



(12) 



Remarkably, this expression does not depend on the 
anisotropy parameter 7 revealing the universal charac- 
ter of the logarithmic negativity at the critical point 
of the LMG model. The fact that a measure of quan- 
tum correlations between macroscopic groups of par- 
ticles is finite at a quantum critical point is also en- 
tirely nontrivial. For example, the mutual information 
1(1, 3) = £ (3, 1 U 2) + £(1, 2 U 3) - £ (2, 1 U 3) diverges as, 
indeed, for the simple case of an equal tripartition, one 
has£(3,lU2) =£(1,2U3) = £(2,lU3). Thus mutual in- 
formation diverges as the entropy, that is, as ^ In [24], 
while it was found to be finite in ID and for short-ranged 
interactions [25, 26]. As I measures all correlations, one 
may conclude that it is the classical part of correlations 
that is responsible for the divergence, whereas the quan- 
tum part (as measured by C) remains finite. We depict 
the behavior of C for such a tripartition and compare it 
with data from exact diagonalization in Fig. 1. 



IV. C IN A BIPARTITION 

Most importantly, when n = r and t 3 = 1 — r, which 
corresponds to the limiting bipartite case r 2 = 0, C di- 
verges at the critical point. Such a behavior is in agree- 
ment with the fact that, for a bipartite pure state, C is 
lower bounded by £ [3] which is divergent at h = 1. To 
analyze this divergence, one expands C in the vicinity of 
the critical point by imposing this bipartition condition 
from the beginning and one obtains : 



£(t, 1 - r,7,/i) 



-iln|^-l| + iln(l- 7 ) (13) 
+ ilnr(l-r) + l + 0(|/ l -l| 1 / 2 ). 



It is interesting to note that this singular behavior is 
exactly the same (up to constant terms) as the one ob- 
tained for other entanglement measures computed in this 
model [15]. This correspondence allows us to straight- 
forwardly extract the finite-size behavior at the critical 




FIG. 1: (Color online) Logarithmic negativity £ as a function 
of magnetic field h for an equal tripartition t\ = ri = T3 = 1/3 
and for two different values of the anisotropy parameter 7. In 
the thermodynamic limit , corresponding to the red (gray) 
lines, at the critical point, C is universal (independent of 7). 
Black lines correspond to numerical data for N = 90, 150, 210 
which match the analytical prediction N = 00 increasingly 
well. Note also that £ vanishes for h — y^y where the ground 
state is separable [12]. 



point by using the same line of reasoning. Indeed, the 
scaling argument introduced in Refs. [11, 12] yields 

C(r, l-r, 7 , 1) ~ \ lniV+1 ln(l- 7 )+± Inr(l-r). (14) 
6 6 2 

In order to check this behavior, we perform exact diag- 
onalization for increasing system sizes at h = 1. As can 
be seen in Fig. 2, numerical data perfectly match the an- 
alytical predictions of the thcrmodynamical limit. Note 
that, in the broken phase, there is an offset of In 2 which 
is due to the fact that the ground state is twofold degen- 
erate. As can be easily understood in the limit h, 7 — > 
for which the finite-size numerical ground state is given 
by a cat-state, this offset is only present for a bipartition 
but does not occur for a tripartition, as can be seen in 
Fig. 1. The expression (14) allows us to add one more 
equivalence of critical scaling laws in the LMG model 
since, at the critical point, we have now 



g „ s ~ £ ~ C, 



(15) 



where Q is the geometric entanglement, S the single-copy 
entanglement, £ the entanglement entropy, and C the 
logarithmic negativity [15]. Of course, it would be very 
valuable to establish the same kind of equivalence in ID 
spin systems for which one already knows that [15] 



N~ 



:£■ 



(16) 



Although a rigorous proof is still missing, our result 
together with recent numerical studies[4, 5] leads us to 
conjecture that in ID critical systems and for a biparti- 
tion, one has £ ^ C. If confirmed, this result may even 
be valid in any dimensions. 
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FIG. 2: (Color online) Logarithmic negativity £ as a function 
of the magnetic field ft for a bipartition t\ = 1/3, tz = 2/3, 
and 7 = 1/2 shown as red (gray) line. In the broken 
phase ft < 1 we plot C + In 2 for a reason detailed in the 
text. Black lines corresponds to exact diagonalization data 
for N = 180, 270, 360. Inset: Scaling of £ with system size 
N = 120, 150 . . . , 1080 from exact diagonalization (black cir- 
cles) approaching a linear dependence on In AT with slope 1/6 
(dashed line) for large N. 



V. THE ISOTROPIC CASE 

Finally, let us discuss the case 7=1 which is triv- 
ially solved since if (7 = 1) commutes with S 2 and S z 
so that the eigenstates are the (permutation-symmetric) 
Dicke states \S, M). For ft > 1, the ground state is fully 
polarized in the z direction (S = N/2 and M = N/2) 
and, consequently, £ = for any tripartition. For ft < 1, 
the nondcgcncratc ground state is still in the maximum 
spin sector S = N/2 but M decreases with ft [12]. The 
isotropic case is thus in a different universality class as 
compared to 7 ^ 1 and it is interesting to compute £ 
in the limit ft — > 1~ . There, the ground state is given 
by \S, M) = \N/2, N/2 - 1) whose logarithmic negativity 
between groups 1 and 3 reads 

£(ri,r 3 ,l,h-> I") = In (l-r 2 + ^ + 4^3). (17) 

This expression strongly differs from Eq. (12) (different 
universality class) while agreeing with the universal char- 



acter of £ in the LMG model at criticality. 



VI. DISCUSSION 

The present study reveals two main properties of the 
logarithmic negativity at a critical point : (i) for a tri- 
partition £ is universal and finite ; (ii) for a bipartition 
£ is universal and diverges as £. To check the general- 
ity of these results, we computed £ in the Dicke model 
[6] for which the ground-state entropy has been already 
computed [24, 27]. This model describes a set of A?" spins 
1/2 interacting with a single-mode bosonic field via the 
Hamiltonian . Thus, if one divides the spins in two parts, 
one can consider two different negativities (spin-spin or 
field-spin). We computed both quantities and we found 
that, contrary to the LMG model, for a tripartition, £ is 
not universal (but still finite) at the critical point. How- 
ever, in the bipartition limit, we found that £ behaves 
also as £ at the transition point. 

These complementary studies of the Dicke and LMG 
models together with ID spin chain analysis [4, 5] lead 
us to conjecture that for a tripartition £ is finite but not 
universal, even at the critical point. Furthermore, in the 
bipartition limiting case, £ and £ behave similarly at the 
transition point. 

A very challenging question would be to check the 
veracity of this conjecture in other spin systems, in par- 
ticular in ID where conformal field theory approaches 
may allow for exact results. 
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